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NORMALIZACION DE LA ESTADISTICA T-MULTIDIMENSIONAL 


por Carlos E. Dieulefait * 


RESUMEN 


Recurriendo a la expresión de la regresión múltiple mediante funciones ortonormales construí- 
das con las variables explicativas, se obtiene la normalización de la estadística t-multidimensional 
descubierta por Ch. W. Dunnett y E. Sobel. Se lograrán así grandes simplificaciones para las tabula- 
ciones correspondientes. Se resuelve también el problema de la estimación de los coeficientes de la 
regresión habitual sin que sea preciso recurrir a la inversión de determinantes. 


ABSTRACT 


On donne la regression multiple exprimée avec des fonctions orthonormales des variables 
explicatives et on presente la normalisation de la formule de la statistique t-multidimensionel decou- 
verte par Ch. W. Dunnett et E. Sobel, ajoutant un nouveau procedé pour la determination de la 
regression multiple habituel sans recurrir aux calculs des determinants. 


1. Sean: 
Yi X]¡ 407 X2¡ +7... FQp Xk¡ +6 =1 + 


Yi =B1 Fi +82 F2¡+... + Bk Fx¡+ e =T;¡ + €j 


dos expresiones equivalentes referentes a un mismo grupo de experiencias. Las y ¡ son los resultados 
de experiencias en la i-ésima prueba o replicación, con i=1,2,...,n>k. Las xy con v= 1,2,..,k 
son las variables explicativas que adquieren las medidas x,¿ en la iésima prueba. Las e, son valores in- 
dependientes de una variable aleatoria normal con esperanzas matemáticas: 


E (e;¡) =0 
E (e; ej) =5;¡0? o 
Las Fy¡ son ortonormales es decir que: 


e o val DK 
1é= 


que se construyen con los valores xy; en la forma: 


Ey; =0y1x]i¡ + 0v2X2i dl O 


Comenzaremos con la primera forma de la regresión. Por la naturaleza de las e¡ se tendrá que 
la probabilidad de tener conjuntamente y, Y2,. . Yn está dada por: 


* Universidad Nacional de Rosario. 


4 ANALES DE LA SOCIEDAD CIENTIFICA ARGENTINA 


n 


Jn 1 
o o 


]= 


El máximo de esta expresión corresponde al mínimo de: 
n 
A 
i=1 
lo que implica que habrá que buscar las soluciones del sistema: 
n 
O = Ti) xy¡=0 Vi nlzo ie 
1= 


Sean a,,2,,.., ax las soluciones o sea los estimadores máximo verosímiles de ,,0%,,.. ,qk. 
Si indicamos con C;¡ a la correspondiente regresión de cálculo, se tendrá: 


C¡=a1X1¡ + 2a2X2¡ +: akXki 


Luego estarán efectivamente realizadas las igualdades: 


n 
(Q) E Ci). =0 =D 
]= 
Haciendo entonces, como es corriente: Y ==; CE) +H(G> E) 
n ó 
y MZ a ES se tiene: 
n= 
(Sa 1 n n 
P=R exp A) explo. 2 (a, a Mae a) Xvi Xji) dy dy2 ...dyn 
20 20% vo = 11 i=1 


con R = (1/0/27 y 


Los n resultados y,, y2, . . , yn tienen n grados de libertad. Por la (2) s* soporta k restriccio- 
nes siendo entonces sus grados de libertad n-k correspondiendo los k g.d.l. restantes a los coeficien- 
tear as. aL. : 

Se tiene así: 


(21 —k — ada 


: 1 An 
exp d a ó 2 (ay —0y) (aj — ay > xviXji) ds? da, da, ... dax 
Mn 1 == 4l 


NORMALIZACION DE LA ESTADISTICA T-MULTIDIMENSIONAL 5) 
con 
(ag 102/4 


r( — mk /2 (n —k)0/2 


siendo A el determinante cuyo elemento (v, j) es: 


n 
2 Xvi Xji 


Haciendo (ay <> ay)/s = ty, sacando a (n-k) s* como factor común en los exponentes e inte- 
grando con respecto as” de cero a infinito se encuentra: 


n n 
2 tyti 2 XyiXj 
A 
Prob (t,,t2,...,tk)=B O O dt,dt, ... dtx 
n — 


vy AT (n/2) Se 


r( + ] A? 


— n/2 


con B = 


Esta fórmula descubierta por Dunnett y Sobel en 1954, está potencialmente llamada en la de- 
cisión múltiple, para dar la distribución de probabilidad de los alejamientos de los coeficientes de la 
regresión de cálculo respecto de los coeficientes paramétricos. 

Uno de los objetivos de este trabajo es el de normalizar la fórmula anterior con las ventajas. 
que además de las teóricas, entendemos se derivarían al emprender un programa amplio de tabula- 
ción. 

2. Trabajando con la regresión ortonormalizada y siguiendo el mismo mecanismo expuesto en el 
párrafo anterior, los estimadores máximo verosímiles by de los f, estarán dados por la fórmula del 
tipo de Euler-Fourier: 


teniéndose así la correspondiente regresión de cálculo: 
C; == b; Ei Al b, E) ¡ + A el biF;¡ 


Pero ahora en lugar del determinante: 


D 
Il 
IMs 


XyiXj¡l| 
Us 


se tendrá: 
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Por ser éste un determinante con los elementos de su diagonal iguales a la unidad y todos los 
otros nulos. 
De este modo resulta la distribución normalizada: 


k 
A 2 7 n/2 
(2) ) e 
Prob (Li, Za... 2) 2 M1 da da da 
mo KN 2 k n=k 
r nm (n=1):? 
2 
Debemos dar ahora la correspondencia entre un dominio D de last,,..., tk a un dominio D de 
las t,, . ., tx. En las aplicaciones estadísticas D y D serán hiperparalelepípedos determinados por las 


dimensiones de sus respectivas aristas. Se deberá tener así: 
Sp Prob. (t, A tx) dt, dt, dtk = $ Prob (Zo Sto 2k ) d2, d2, Sn Uk 
Debido a que: 


(3) a1X + 2a2X2¡ +7 + agxx¡ = dy Fj¡ +b2F)¡+>::+0x Fri 


ax +mx +: Tax PPPoE 


se deduce que: 


xi tt Aa FB Be. 2 E 


Luego teniendo en cuenta la (1) resultan las siguientes fórmulas de pasajes: 


ty = 2y Oyvy ARE O: +2 Uk y 


y las relaciones simétricas: 


Ze ty Y “e tv+i Kvev+1 + tk Av. .k 


donde: n 
Aj¡= 2 On Xi 


1= 


3. La fórmula (3), teniendo en cuenta la (1) es la que nos permite determinar los coeficientes de 
la regresión habitual en función de la regresión de cálculo normalizada. 
Se tiene, de este modo: 


dv =DyO0yv 1 Dy+10u+r.v to... +? OKOk-y 


con lo que se puede prescindir por completo de todo cálculo de inversión de determinantes como 
es corriente. 


Referencia. 


Ch. W. Dunnett y E. Sobel. Rev. Biometrika 1954: A bivariate Generalisation of Student's t-distribution, with ta- 
bles for certain special cases. : 
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ESTABILIDAD DE SOLUCIONES DE 
HBr - ESCOPOLAMINA DE USO FARMACEUTICO 


por Juana HUGHES! y Patrick MOYNA? 


(1) Dexin S.R,L., Colonia 1181, Esc. 305, Montevideo, Uruguay. 


(2) Cátedra de Farmacognosia, Facultad de Química, Montevideo, Uruguay 


RESUMEN 


Se realizaron ensayos acelerados de descomposición de una solución de escopolamina en buffer fosfato 
pH 4,0. La concentración se determinó por cromatografía gaseosa. Se calculó la k de la velocidad de reacción, 
y su variación con la temperatura. A partir de la ecuación de Arrhenius se estimó el tiempo en que a 20 *C la con- 
centración se reducirá al 90% de su valor teórico. 


ABSTRACT 


Stability of scopolamina in buffer solution pH 4,0 was investigated, k for the rate of reaction was determi- 
ned, together with its relationship with temperature, From the Arrehenius equation the time taken to lower the 
concentration to 90% of the initial value at 20 “C was calculated. 


INTRODUCCION 


La escopolamina se usa como anticolinérgico desde hace más de un siglo. Una de las formas 
farmaceuticas en que se prepara es en inyectable (soluciónacuosade bromhidrato de escopolamina) (? ). 
Al igual que otros alcaloides del grupo del tropano, es un éster del ácido trópico (?). Se hidroliza en 
solución acuosa a velocidades que dependen fundamentalmente del pH y la temperatura, y que pue- 
den ser muy rápidas (?) y (*). La cinética de esta reacción ha sido poco estudiada. Para ello, es 
necesario emplear una técnica de dosificación que sea específica para la escopolamina en presencia 
de sus productos de degradación. Si bien hay descritos varios métodos para dosificar la escopola-- 
mina (5), ($), y ( Y *£. ahí incluídas) muy pocas sirven para estudiar la cinética de su hidrólisis. La 
técnica introducida por la U.S.P. XVHI por G.L.C. es la más adecuada por su sensibilidad, espe- 
cificidad y porque permite hacer una serie de medidas suficientemente grande en plazos cortos. 


MATERIALES Y METODOS 


—Solución estandard (SE) 


Se preparó una solución de HBr — Escopolamina estandard (5 mg/ml) en buffer de ácido 
cítrico — Na, HPO, 0,02 M pH 4,0. 


—Procedimiento 


La solución estandard se fraccionó en ampollas de 1 ml (vidrio USP XVIII tipo I) y se dividió 
en grupos que se termostatizaron a 40%, 80” y 108*C. Al completarse el tratamiento térmico, las 
amnpollas se dejaron enfriar. De cada una se tomaron 0,8 ml de solución que se llevó a pH 9,0 con 
0,8 ml de solución buffer NaH2PO¿ — NaOH 0,1 M pH 9.0. Se extrajo con 1 ml de CHClz que 
contenía 0,2% de miristato de metio (estandard interno). Se repitió el proceso de extración con 0,8 
ml de la solución estandard sin tratamiento térmico. 
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—Cromatografía gaseosa 


Se inyectó 1 ul de extracto clorofórmico de SE y solución problema sucesivamente en un cro- 
matógrafo de gases. Se utilizó un G.L.C. PYE Unicam 104 con detector FID dual, integrador PYE 
Unicam DP 80 y registrador PHILLIPS DM 8100. Se empleó una columna de vidrio (largo 0,90 m, 
diámetro interno 4 mm) con SE 30 al 5% (sobre Diatomite CAW DMCS 100 — 120 mallas) como 
fase estacionaria. Se trabajó con 200 *C de temperatura del detector. Flujo de N, : 30 ml/min. 


—Cálculos 


La concentración de escopolamina sin descomponer fue determinada para cada valor de tiem- 
po y temperatura, tomando 3 ampollas diferentes, dosificando 3 veces cada una de ellas. Se calcula 
la concentración como: 

Cs. concentración de SE 
o hms-. altura del miristato de metilo en cromatograma de SE 
O ETE hp. altura de Escopolamina en cromatograma problema 


hes X h 
ved hs. altura de Escopolamina en cromatograma de SE 


hip: altura de Miristato de metilo en cromatograma problema 


Las pendientes e intercepciones de las rectas se obtuvieron en una minicalculadora Texas Ins- 
trument SR 5] A. 


RESULTADOS Y DISCUSION 


El pH de la solución se seleccionó como el más conveniente para una solución de escopolamina 
teniendo en cuenta su influencia en la velocidad de descomposición (?) y (*). 


Para cada temperatura se calcularon las ecuaciones de las rectas obtenidas al graficar C vs t, 1n 
C vs t y1/C vst(?). En todas ellas se encontró que los resultados mejor alineados eran los de ln Cy. t 
o sea la ecuación de primer orden. 

Luego se calculó en forma analítica la variación de ln k vs 1/T para obtener los parámetros de 
la Ec. de Arrhenius: 


A H 
In k = + In A 
RT 


TABLA I 


| a a EA ep 


OOO Al O 
A see 1 c4 
13 0380) 08 
A O 


43 82,5 | OS 
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- TABLA UH 


Tp CO) | Haas) 


| 708 


GRAFICA 1 


Ink 


> 4 (210?) 
VE 1 1 ¡ +] 
20) u 


2.0 2.5 3.0 


Considerando que la temperatura media sea de aproximadamente 20 *C se puede calcular que 
el tiempo en el cual la concentración inicial de escopolamina en solución acuosa a pH 4,0 se reduce 
al 90% del valor teórico es de 2,5 años. 
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ESTUDIOS POR SIMULACION ANALOGICA ELECTRONICA 
DE LA POBLACION DE UNA ESPECIE ICTICOLA EN UN LAGO 


Por C. C. SOTERAS?, C. A. LEGUIZAMON!, V. M. BRICEL]J?. 
J. L. CHURRUARIN?, A. C. REPETO?, S. TOKUSHIRO? 


ABSTRACT 


The problem is analized for populations of a fish of 2, 3, ..., n years. Baside all the ecologycal facts that 
normally have an influence in the number variation, it is also considered the effect of the setting of an Nuclear 
Power Plant, which uses big quantities of water for cooling. This last thing has two important ways of acting over 
the fish: thermic shock and mechanical damages. 

The object of this study is the consideration of these alternatives in order to determine the magnitude and 
moment of influence. As it is necessary to keep the number of fishes never under the optimum number it is possible 
to recognize what action to follow over the ecologycal system considering these situations. 

Whith the electronic analog simulation it can be constantly controlled the variation in the number of 
fishes. Thus, after an initial knowledge of the need values for the different equations that are used, it will only be 
necessary to make fishes census from time to time, each time more separated, in order to regulate the eventual 
deviations of the system. 


RESUMEN 


El problema es analizado para poblaciones de un pez de 2, 3,...,n años. Además, de todos los factores 
ecológicos que normalmente intervienen en la variación del número, son considerados también los efectos produ- 
cidos por el emplazamiento de una Central Nuclear, la cual utiliza elevados caudales de agua para enfriamiento. 
Esto último, provoca dos importantes acciones sobre el pez: shock térmico y daños de origen mecánico. 

El estudio se basa en la consideración de estas alternativas, determinando la magnitud y momento de 
incidencia de las mismas y dado que es necesario mantener un número de peces que no debe ser inferior al óptimo, 
se puede reconocer la acción a seguir sobre el sistema ecológico. 

La simulación analógica electrónica permite contar con una metodología óptima para el control permanente 
de la variación del número de peces. Se prevee que después de un conocimiento inicial de los valores para imple- 
mentar las ecuaciones a usar, sólo será necesario censos de peces cada vez más temporalmente espaciados para 
ajustar las eventuales desviaciones del sistema. : 


INTRODUCCION 


Este trabajo pretende introducir la tecnología de la simulación análogica electrónica para 
efectuar el control permanente de la variación de las poblaciones de una especie de un pez en un 
lago. El emplazamiento de una Central Nuclear a orillas de un lago utiliza grandes caudales de agua 
para refrigeración, pudiéndose producir una disminución en el número de peces (1) tanto por 
elevación de temperatura del agua (2, 3,4, 5) (shock térmico) como por daños de origen mecánico. 

Los peces son organismos poiquilotermos y la temperatura juega un papel muy importante 
como regulador de su fisiología, por controlar el ciclo reproductivo, la velocidad de digestión, la 


NOTAS 


(1)  Biomatemática, Gerencia de Investigaciones, Comisión Nacional de Energía Atómica y Depar- 
tamento de Matemática, Facultad de Ciencias Exactas y Naturales, U.B.A, 


(2) Comisión Nacional de Energía Atómica. 


(3)  Biomatemática, Facultad de Ciencias Exactas y Naturales. U.B.A. 
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tasa de respiración y la incidencia de las enfermedades. La temperatura no afecta a los organismos 
tanto por su valor absoluto como por el salto térmico. Pequeños saltos de hasta 0.5” C, pueden ser 
detectados por los peces. Los efectos del incremento de temperatura sobre la población ictícola se 
pueden dividir en directos e indirectos. Los indirectos se observan principalmente sobre la fisiología 
y el metabolismo produciendo cese de la reproducción, afectando a las fuentes de alimento, 
favoreciendo el crecimiento y reproducción de los parásitos y teniendo como consecuencia un 
aumento en la tasa de infección. La muerte por acción directa se observa raras veces en la naturaleza 
ocasionando más efectos en los estadíos más tempranos, los cuales son más vulnerables. como así 
también durante la incubación produciendo malformaciones o un cese dei desarrollo. 

La metodología es impuesta por criterios concretos de análisis; esto es, controlar las pobla- 
ciones de distintas edades partiendo de datos iniciales obtenidos, y cuyas desviaciones se ajustarían 
per una sencilla técnica experimental como es el censo periódico de peces. 

Estos datos iniciales deben ser obtenidos experimentalmente para implementar la ecuación 
diferencial adoptada, cuya solución bajo distintas condiciones, se pretende hallar. 

Sobre la base de un número mínimo de peces de cada población, el análisis permitirá determi- 
nar momentos y cantidades de nuevos peces a sembrar en el lago, así como, introducir los datos 
emergentes de este tipo de estudios ecológicos, en el conjunto de otros criterios para la regulación 


del funcionamiento de la Central. 
ANALISIS DEL PROBLEMA 


La figura 1 es ilustrativa del problema. La Central Nuclear se puede encontrar emplazada 
de tal forma que el agua de refrigeración que utilice tome un valor de temperatura Ts a la salida 
de los intercambiadores de calor, la cual después de recorrer el canal de drenaje ingresa nueva- 
mente al lago a temperatura Tc. enfriándose naturalmente hasta la temperatura T en superficie 
y Tf a nivel de la toma de agua. Al solo efecto ilustrativo, tenemos: 


AT =Tf-— Ts<7*C 
ISTMO = 00516 


FIGURA 41. El lago y el emplazamiento de la Central Nuclear. A. Central Nuclear; B. Toma de agua a aproximadamente 15 m de 


profundidad. E. Canel de drenaje del agua que retorna al lago. T. Temperatura én la superficie dei agua del lago. 


T;. Temperatura del agua en la toma. T¿= Temperatura del agua a la salida de la Central . (Ts = Y + A1;At= salto térmico). 


T,= Temperatura del agwa a la solida del canal. (TFT, - 0,5%). 


lo cual depende principalmente de la región de emplazamiento, características técnicas de la Cen- 
tral, profundidad de la toma de agua y de la época del año. 

La profundidad de la toma de agua es importante, no sólo por cuanto al ser mayor, menor es 
Tf lo cual favorece la capacidad natural de refrigeración del lago (a AT cte, Tc es menor), sino que 
también, dado que el pez no habita en regiones carentes de oxígeno, entonces es posible un efecto 
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nulo sobre éste, si la toma es efectuada a gran profundidad, siempre que además, el pez no sea arras- 
trado por la succión. A este respecto, es necesario hacer notar que provista la boca de malla protec- 
tora o equipos electrónicos de repulsión del pez, puede ser nulo el ingreso de éste; o bien, que sólo 
ingresen peces de pequeño tamaño. 

Las hembras desovan en lugares del lago de poca profundidad, durante dos ocasiones en el 
año, con una cantidad de huevos dependiendo principalmente de la edad de la hembra y la época 
de desove. El número de huevos sembrados pueden ser aleatoriamente disminuídos, tal como por 
simples variaciones en el nivel del lago, posible daño mecánico de la Central por efectos secunda- 
rios, etc., entre otros. Estos importantes motivos, sumándose a otros inherentes a la compeljidad 


ecológica: factores físicos (luz, temperatura, oxígeno), químicos (concentración de nutrientes en 
el lago y desde ríos y arrollos) y orgánicos (fitoplancton, zooplancton) se interaccionan con las 


conocidas sucesivas etapas o estadíos del crecimiento biológico de la especie: 

Alevinos > Juveniles > Adultos 

Razones prácticas del análisis, tal como el hecho dado por el interés turístico sobre poblacio- 
nes de edad de y mayor de dos años y también por la acción de los efectos de la Central según la 
edad de cada población, hacen que el estudio se efectúe considerando la variación en el número 
de peces de cada edad y en el período de un año. De hecho, interesará mantener un número no 
inferior a un valor M¡(i, 2, 3,.. . ,n) de peces de edad i al final de cada año, a los efectos de 
preservar la especie en el lago. 

La edad del pez es determinada por la interpretación y lectura de marcas de crecimiento en 
estructuras óseas, tales cómo escamas u otolitos. La escama característica de un pez, que se ilustra 
en la figura 2, es una estructura que se forma por oposición de delgadas placas óseas, y que comien- 
za a generarse recién al finalizar el consumo del saco vitelino de la larva. En la escama perteneciente 
a un pez de regiones templadas, que sufre variaciones estacionales en su ritmo de crecimiento, es 
posible observar cada año una banda de líneas de crecimiento muy próximas (que corresponden a la 
estación fría, de crecimiento nulo o retardado), seguida por una banda de líneas espaciadas (que 
corresponde al período de crecimiento rápido de los meses más cálidos del año). Así el recuento 
de bandas del primer tipo, que se denominan anillos de crecimiento o annuli, permite estimar la 
edad del pez y correlacionarla con la talla. Si se supone el caso de una especie que posee dos perío- 
dos de desove al año, por ejemplo primavera y verano, un observador experimentado podría diferen- 
ciar los peces provenientes de cada camada, por la distancia entre el foco de la escama y la primer 
marca invernal. sin embargo, a los fines de una mayor significación, para el modelo se consideró 
como si los peces ton una sola marca invernal fueran todos uniformemente de una dada edad, inde- 
pendientemente de la distancia de dicha marca al foco. 


Campo anterior 


Campo lateral 


SN: Us Anillo de crecimiento 
2 yl 
y) 
Y 
== $ 1 Radios 


ES), 


Foco 


Campo posterior 


FIGURA 2. Esquema de la escama de un pez con dos anillos de crecimiento 
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De esta forma, se analizan las poblaciones desde el 1? de setiembre (con la edad recientemente 
determinada en sus escamas) hasta el 31 de agosto del año siguiente. 

Considerados los dos períodos de desove con medias en marzo y octubre de cada año, es inme- 
diato que una población de edad 2 años al 1% de setiembre del año m, corresponde a huevos deso- 
vados en octubre del año m-2 y en marzo del año m-1. 


Para edades 1,2,3,...,n años, planteamos la ecuación diferencial 
dXi E A 
— =XoW + ke UV Xi —KU xi (1) 
dt 
donde: 


Xo() = número inicial de peces de la población i. 
Ke" = proporción de peces de edad i que ingresan por unidad de tiempo desde ríos y arroyos. 
kK0U = proporción de peces de edad i que se eliminan por unidad de tiempo. 


El valor de K“? es dado por: 


k0-"5 10 (2) 
j= 1 
donde cada k,' es dado por la nomenclatura. 

Unas pocas palabras deben ser dichas de algunos de estos coeficientes. Tanto kel como k; | que 
determinan la proporción de peces que ingresan y egresan del lago, son tomados como valores cons- 
tantes, aunque son variables según la época del año (por ejemplo, el período en el cual las hembras 
buscan el desove fuera del lago). 

El valor de k»! sin ser constante (aunque así es tomado) no es significativo. Los valores de K3! 
y ky¿! requieren una compleja experimentación, aunque kx! es un valor escasamente significativo; 
mientras que ks! si bien es tomado constante, es necesariamente variable en el tiempo, en función de. 
la época de pesca deportiva o por la época de pesca furtiva, o semicomercial, que también podría existir. 

Cada valor de kg determinado por la acción mecánica de la Central es muy dependiente de la 
edad del pez ya que las edades mayores se supone que no ingresan a través de la malla que debería 
existir en la boca de toma. También k¿! depende del régimen de la Central. El valor de k, 1 de peces 
muertos naturalmente tiene un valor escasamente significativo y es más importante en poblaciones de 
mayor edad. Este coeficiente k,1en realidad debe presentarse ligado a kg i que considera la variación de 
la temperatura del agua en el lago. Considerando una variación de tipo senoidal en la temperatura T 
(t) promedio en el lago (las características de la sinusoide es variable según los años), para kg1se toma 


kg = kg (T(t) — To) 3) 


donde To =* 22 C, temperatura umbral de supervivencia del pez. Así, para T (t) = To, el efecto de 
muerte es dado por k+. 

Respecto del shock térmico, debe consignarse que los efectos térmicos de este tipo no son sólo 
dependientes de la especie, sino fundamentalmente también del valor AT del shock y de la tempera- 
tura T a la cual se encuentra normalmente el pez. Dado que el presente trabajo tiene otros fines dis- 
tintos a la determinación precisa de cada factor incidente, se toma para ko 


ko Y =k, 0 (T (t) + AT To! (4) 
donde AT = Tf — Ts; To =22"C. 
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ASPECTOS DE LA SIMULACION Y RESULTADOS 


A los efectos de mostrar menores detalles de la simulación del problema, es útil analizar el dia- 
grama simbólico general ilustrado por la figura 3. Este diagrama con algunas pocas modificaciones de 
crtierio práctico, fue el utilizado para obtener los datos y curvas a las que posteriormente nos refe- 


riremos. 


Sobre la base de las consideraciones anteriores, se tiene la siguiente ecuación (5), básica para la 


simulación 


FIGURA 3. Diagrama simbólico general para simular las variaciones de poblaciones de distinta edad de un pez 


en un lago. 


dXi o AS E (k, ae k, 04 k, Oy KO + k, 1) Ma — (Ka 4 sE k¿)) xo 


ka D (T (t) - To) Xi—k9% (T (t) + AT —ToY Xi (5) 


Refiriéndonos a la figura 3, la ecuación (5) es implementada de la siguiente forma 


La variable Xi se resuelve mediante el integrador 1. 


El valor inicial de peces Xo! es dado por la condición de entrada (IC) multiplicada por el 
factor del potenciómetro respectivo. 


El factor a, del potenciómetro 2 tiene en cuenta al valor ke! debido al ingreso de peces de 
edad i desde ríos y arroyos. 


Toda posible combinación de valores expresados por los coeficientes del tercer término del 
segundo miembro de (5) es dado por «3, aunque k», (1), k3(*) y k, (1) pueden ser nulos o esca- 
samente significativos según la edad del pez. La facilidad de operar e imponer valores sobre 
los potenciómetros, permite tener en cuenta, muy fácilmente, las anteriores alternativas. 

El efecto de muerte del pez por pasaje a través de la zona caliente del lago (k4”) y el debido 
a la disminución de peces por el pasaje del mismo a través del circuito de agua de refrigera- 
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ción (ké) por ser propios del funcionamiento de la Central, se incluyen bajo el valor del 
factor 4 del potenciómetro 4. 

— El efecto de muerte por temperatura promedio del lago se obtiene a través de la línea del 
diodo D1. La función senoidal característica de la variación de la temperatura del agua del 
lago durante el año, se compara a través del pasaje o no de corriente por D1, con la fuente 
El cuyo valor representa la temperatura To = 22 2C que es la umbral de supervivencia del 
pez. El valor positivo obtenido cuanto T (t) > To, se multiplica con Xi en el multiplicador 
1, para luego Q% tener en cuenta ka). 


— El shock térmico es simulado por la línea del diodo D2. La temperatura del lago variable 
senoidalmente durante el año, se le suma un valor dado por la línea del potenciómetro 7, el 
cual tiene en cuenta el salto término producido por la Central (AT =Tf— Ts). Aquí es intere- 
sante notar la igualdad -o no de T (t) con Tf que puede hacerse, según las características y 
profundidad de la toma de agua. El diodo D2 permite el pasaje de corriente cuando E2 
representado To = 22” C es superado. El multiplicador 2 eleva al cuadrado y él multipli- 
cador 3 incluye el valor de Xi. El potenciómetro 8 tiene en cuenta e La salida ingresa al 
integrador 1. 


Sobre las bases generales del circuito correspondiente al diagrama simbólico de la figura 3, se 
analizan algunos de los casos obtenidos, a través de las tablas y figuras consignadas. 


TABLA 1 


FIGURAS 
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TABLA 2 


FIGURAS 
6-a 6-b 6-c 6-d 7-a 7-b 7-cC 7-d 
X, 1.000 1.000 1,000 1.000 1,000 1,000 1,000 1,000 
6 RE 2336 2396 DE 2300 23,76 390 
E 50 230 506 50 5 E sc 
AT 156 ae 720 za E 76 7*C 
a A 0,021 0,021 0,021 0,200 0,200 0,200 0,200 
00917 0,017 0,017 0,017 0,090 0,090 0,090 0,090 
01? 0,012 0,012 0,012 0,060 0,060 0,060 0,060 
O 0,010 0,010 0,010 0,030 0,030 0,030 0,030 
o 012 0,012 0,012 0,012 0,060 0,060 0,060 0,060 
Eo ole 0,016 0,016 0,016 0,060 0,060 0,060 0,060 
a 0.067 0,047 0,067 0,067 0,300 0,300 0,300 0,300 
k 0 ss pa 0,002 0,002 El E 0,002 0,002 
k 0 se E 0,008 0,142 pa s 0,008 0,142 
ES = E 0,010 0,144 de ES 0,010 0,144 
O 0/0 1a 0,012 0,012 0,012 0,012 0,012 0,012 0,012 
00 0,001 0,001 0,001 0,001 0,001 0,001 0,001 
TABLA 3 
FIGURAS 
st... 8 9 10 
X, 1.000 1.000 1,000 
aa 2380 25% 25*C 
e 756 7ec 7% 
AT E E E 
5 0,005 0,005 0,005 
k 0 0,017 0,017 0,017 
k 0 0,012 0.012 0,012 
k, 0 0,010 0,010 0,010 
k.0 0,012 0,012 0,012 
k, 0 0,016 0,016 0,016 
0% 0,067 00,067 0,067 
x 0 E sl A 
O - - - 
k, 6 0,015 0,020 0.025 
k 0 0,005 0,010 0,015 
COMENTARIOS 


El sistema ecológico como el planteado, o sea la variación del número de individuos en pobla- 
ciones de distinta edad de un pez en un lago, admite la consideración de mayor número de factores, 
más complejidad en las interrelaciones que aquellas aquí consideradas, y además, análisis más pro- 
fundos de muchos coeficientes de la ecuación (5) que no son simplemente constantes. Ellos pueden 
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ser variables en el tiempo (muy en especial ke! an y key o bien variables en función de otras varia- 
bles dependientes del mismo sistema ecológico. A estas consideraciones, deberíamos adicionar otras 
de origen tecnológico, tal que desde ambos criterios, no podría el problema, llevado a la práctica, 
ignorarlos. No obstante lo anterior, la metodología pretende afirmarse sobre la practicidad de las 
determinaciones, como por ejemplo, la corrección periódica de las desviaciones de las curvas por 
medio de censos. Esto evidentemente nace de una viabilidad de este tipo de estudios, limitados por 
factores económicos. 

Respecto de la implementación de las soluciones, las figuras e leccin la mayor disminu- 
ción de peces por: el efecto de la temperatura en época estival, el efecto debido al shock térmico, 
(también más significativo en el verano), y los efectos mecánicos del pasaje del pez por la Central. 
Estas mayores. 0, menores variaciones según el período anual, permiten determinar el momento 
de la siembra de nuevos peces. La cantidad de esta última, será dada por el número mínimo sobre 
cada población a mantener. De hecho, se cumplirá la igualdad Mi = Xo y esto permite el estudio 
de un pez, en distintas edades bajo los mismos efectos. Esto, correspondería a un análisis preventivo. 
De aquí surge la segunda alternativa, cual es tratar los distintos valores iniciales Xo! = 2, 3, ,n) 
existentes en un lago en setiembre de un dado año, y valorar su variación según todos los efectos 
naturales y o: que afectan, desigualmente, a todas las poblaciones. 


NOMENCLATURA 


Xo!: Cantidad inicial de peces de edad i. 

kel: Proporción de peces de edad i que ingresan por unidad de tiempo desde ríos y arroyos. 
k: Proporción de peces de edad i que emigran del lago por unidad de tiempo. 

k!: Proporción de peces de edad i predados por unidad de tiempo. 

ka: Proporción de peces de edad i, muertos por desnutrición por unidad de tiempo. 


ka: Proporción de peces de edad i, muertos por pasaje a través del gradiente de temperatura 
que se forma en el lago, a la salida del canal de drenaje, por unidad de tiempo. 


ko: Proporción de peces de edad i, extraídos por pesca por unidad de tiempo. 
ké: Proporción de peces de edad i, muertos mecánicamente por la Central. 
ko: Proporción de peces de edad i, muertos naturalmente por unidad de tiempo. 


kg: Proporción de peces de edad i, muertos por variación de la temperatura promedio en el 
lago por unidad de tiempo. 


ko!: Proporción de peces de edad i, muertos por el shock térmico. 
10% Tiempo. 

To: Temperatura umbral. 

AT: salto de temperatura. 

Tmax: temperatura anual máxima en el lago. 


Tmín: temperatura anual mínima en el lago. 
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FIGURA 6. Condiciones dedes por la iabla 2 
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FIGURA 7. Condiciones dedas por la tabla 2 
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FIGURA 8. Condiciones dades por la tabla 3 


(a) Función senoidal de amplitud IArepresentativa de la vorioción anual de la temperatura del lago) 
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FIGURA 9. Condiciones dadas por la tabla 3 


(a) Función senoidal de amplitud 92 (representativa de la variación anual de la temperatura del lago) 
(b ) Número de muertos por choque térmico 


(c ) Número de muertos por temperatura del lago 
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FIGURA 10 Condiciones dadas por la tabla 3 


(a) Función senoidal de amplitud 9? (representativa de la variación anual de la temperatura del lago) 
(b) Número de muertos por choque térmico 


(c) Número de muertos por temperatura del lago 
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ESTUDIO CON MICROSCOPIO ELECTRONICO DE BARRIDO DE LA 
EVOLUCION DE LAS CELULAS CILIADAS ECTODERMICAS EN EMBRIONES 
DE BUFO ARENARUM 


Por JORGE HERKOVITS* y ARMANDO PISANO** 
RESUMEN 


Las cilias comienzan a crecer sobre la superficie de algunas células del ectodermo durante el estadio de néu- 
rula (estadios 13-14). El número de células ciliadas se incrementa rápidamente de tal manera que en el estadio 
17-18 muchas células ciliadas contactan entre sí, al mismo tiempo que el número de cilias per célula también 
aumenta. En estadios subsiguientes (estadios 19-23) se observa que las células ciliadas aparecen gradualmente más 
separadas entre si por células glandulares y no se forman nuevas células ciliadas. En los estadios 23-25, las células 
ciliadas presentan gradualmente menor número de cilias y finalmente, en la superficie apical de estas células, que- 
dan pliegues que presentan una imagen reticular. Se cumple así el ciclo de presencia de cilias en el ectodermo de 
Bufo arenarum. ; 


ABSTRACT 


Cilia start to grow from the free surface of some ectodermal cells during neuruia stage (stage 13-14). Then 
the number of cellated cells increases greatly and as a result during early tail bud stage (stages 17-18) many 
ciliated cells are in contact each other; in the other hand the number of cilia per cell also increases. During subse- 
cuent development (stages 19-23) the ciliated epidermal cells appear gradually more separated by secretory ecto- 
dermal cells. By stages 23-25 regression of cilia can be detected and the free surface of these cells exhibits folds 
forming a reticularlike structure. The observations confirm that during amphibian development there is a cycly 
by which cilia appear aná disappear from the ectodermal surface. 


INTRODUCCION 


A partir del estadio de néurula, el epitelio ectodérmico de embriones de anfibios presenta cé- 
lulas ciliadas (Assheton, 1896; Pflugfelder y Schubert, 1965, Billet y Gould, 1971; Billet y Courte- 


nay, 1973; Kessell y col., 1974). Las células ciliadas cumplen un ciclo de presencia en el tejido 
ectodérmico motivo por el cual resulta de particular interés estudiar los cambios morfológicos que 
acontecen en la superficie aplical de estas células durante el desarrollo embrionario. 

Para el análisis de las características de la superficie celular, el Microscopio de Barrido (MEB) 
posibilita una útil combinación de resolución, contraste y profundidad. En efecto, estudios realiza- 
dos siguiendo procesos morfogenéticos que acontecen in vivo e in vitro (Tarin, 1971; Billet y Cour- 
tenay, 1973; Kessel y col., 1974; Ben Shaul y Moscona, 1975; Monroy y col. 1976; Herkovits, 1976 
a; b; Herkovits y col., 1977; etc.) han demostrado que cambios morfológicos que acontecen en la 
periferia celular se encuentran asociados a procesos de migración, interacción y diferenciación ce- 
lular. ¡ 

En el presente trabajo, se estudió con MEB el patrón témporo-espacial de la formación, distri- 
bución y desaparición de las células ciliadas sobre la superficie ectodérmica de embriones de Bufo- 
arenarum. Asimismo se presentarán algunas características de la superficie apical de las células ecto- 
dermicas glandulares. j 


* La dirección actual del Dr. Herkovits es “Instituto de Biología de la Reproducción y Desarrollo Embrionario”. 
Universidad Nacional de Lomas de Zamora cc 95(1832). Lomas de Zamora, Prov. de Buenos Aires. 


** Laboratorio de Investigaciones Embriológicas, Consejo Nacional de Investigaciones Científicas y Técnicas, 
Paraguay 2155 Buenos Aires. Investigador del CONICET. 
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MATERIAL Y METODO 


Los ovocitos de Bufo arenarum fueron obtenidos por inyección intraperitoneal de una suspen- 
sión de hipófisis homologa mantenida según Pisanó (1957). Los cordones ovígeros extraidos del ovi- 
saco y colocados en cápsulas de Petri, fueron fertilizados in vitro con una suspensión de espermato- 
zoides. Al alcanzarse el estadio de blástula se procedió a disolver la cubierta gelatinosa que rodea a 
los embriones, usando ácido tioglicólico al 2% (pH 7,2). Para el estudio del ectodermo con el MEB, 
fueron fijados embriones entre los estadios 13 y 25 (ver tabla normal de Del Conte y Sirling, 1952) 
con glutaraldehido al 2,5% en buffer fosfato, deshidratados en una serie ascendente de etanol y se- 
cados de acuerdo con la técnica del punto crítico. Los especímenes secos fueron recuiertos con oro- 
paladio evaporado y observados con un MEB Jeol JSM-U3 que fue operado entre 5 y 15 kW. 


OBSERVACIONES 


Las cilias aparecen sobre la superficie ectodérmica de embriones de Bufo arenarum entre los 
estadios 13 y 14, es decir entre el estadio de placa neural y la formación del canal neural. Durante 
este período temprano las células ciliadas se encuentran muy separadas entre si por células ectodér- 
micas glandulares (Fig. 1) y cabe mencionar que presentan una superficie apical menor y de forma 
más circular que las células glandulares que son poliédricas. Las células que constituyen el ectoder- 
mo neural no desarrollan cilias. El crecimiento de las cilias pertenecientes a una determinada célula 
no es sincrónico, como tampoco lo es el crecimiento ciliar entre distintas células y/o embriones de 
un mismo estadio; no obstante al alcanzarse el estadio 16 la superficie apical de numerosas células 
está totalmente cubierta por cilias. La formación de nuevas células ciliadas continua hasta el estadio 
17-18 de tal manera que numerosas células de este tipo se encuentran contactando entre si (Fig. 2). 
Cabe mencionar que mientras los esbozos ciliares en algunas células se encuentran en la parte central 
de la superficie apical, en otras se observa que su distribución es más uniforme sobre toda esta super- 
ficie (Fig. 1 y 2). En estadios ulteriores (v. g. 19-23), las células ciliadas se encuentran separadas 
entre sí por células glandulares, de tal manera que como regla general se puede decir que cada célula 
ciliada está rodeada por células glandulares. Además toda la superficie del embrión, con excepción 
de pequeñas áreas (v.g. la superficie de las ventosas, la placoda cristalineana, parcialmente el esto- 
modeo) presentan células ciliadas (Fig. 3). 

La regresión de las células ciliadas se produce a partir del estadio 23-24, siendo primero mani- 
fiesto en las regiones laterales de la cabeza, del cuerpo y de la cola del embrión. La desaparición de 
estas células continúa en forma asincrónica (al igual que su aparición), pudiendose observar en el es- 
tadio 24 en una misma región (v.g. aleta caudal) células que aún conservan la mayor parte de su 
población ciliar cerca de otra célula que presenta solamente unas pocas cilias de longitud variable 
pero siempre mucho más cortas que durante el período de su máximo desarrollo. En otra célula cer- 
cana se observa que las cilias han involucionado casi totalmente (Fig. 6). En las células que presentan 
signos marcados de involución ciliar, se puede observar la presencia de numerosos pliegues que se 
anastomosan entre si dando la imagen de un retículo sobre la superficie apical de la célula. Ulterior- 
mente los pliegues también desaparecen cumpliéndose así un ciclo de existencia de las cilias sobre la 
superficie ectodérmica de embriones de Bufo arenarum. 

Las células glandulares presentan numerosas microvellosidades que en muchos casos se anasto- 
mosan entre sí, observándose también pequeños gióbulos (ca 0,8 um de diámetro) y esporádica- 
mente se encuentran también pequeños orificios en la membrana celular (Figs. 1, 2,3, 4). 

La superficie de contacto apicolateral entre células ciliadas y/o glandulares sobresale sobre la 
superficie epitelial lo que permite visualizar con mucha nitidez los límites celulares (Figs. 1, 2, 3, 4). 
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Examinando cortes ultrafinos con el microscopio electrónico de transmisión se pudo comprobar que 
la región de contacto intercelular a nivel del borde apico lateral presenta numerosos desmosomas por 
debajo del complejo de unión apico lateral (Fig. 5). 


DISCUSION 


El tiempo de aparición y la distribución de las células ciliadas pertenecientes al ectoaermo de 
embriones de Bufo arenarum es similar al observado en otros anfibios (Billet y Courtenay, 1973; 
Kessell y col., 1974) siendo el aspecto general de las cilias similar al observado en otros tejidos (Bar- 
ber y Boyde, 1968; Sleigh, 1962). Con respecto al origen de las células ciliadas se ha sugerido que los 
precursores se encuentran por debajo de la capa ectodérmica superficial y que la formación de cilias 
acontecería cuando esas células se desplazan hacia el estrato superficial (Billet, 1968; Billet y Cour- 
tenay, 1973). Coincidentemente con esta sugerencia hemos observado que en estadio de néurula las 
células que desarrollan cilias presentan una superficie apical mucho menor que las células glandulares 
vecinas, lo que haría suponer su desplazamiento desde la capa inferior. 


Considerando que estructuras como los folículos pilosos (v.g. en la piel de la oveja, Claxton, 
1964), los pelos de la superficie de los insectos (Lawrence, 1970) las estomatas en las hojas (Bunning, 
1965) se encuentran alejados entre sí, siguiendo una distribución espacial específica (Dick Gordon, 
1971), se postuló que debería existir un mecanismo inhibitorio capaz de impedir el desarrollo de ele- 
mentos similares en contacto entre sí (Lawrence, 1970; Wigglesworth, 1940). No obstante y si bien 
como regla general se puede decir que las células ciliadas se encuentran separadas entre sí por células 
glandulares, es de destacar que en los estadios 17-18, cuando la densidad de células ciliadas parece 
máxima, numerosas células de este tipo contactan entre sí. Esta observación permite sugerir que no 
existiría un fenómeno de inhibición “lateral” por parte de las células ciliadas para la formación de 
otras células del mismo tipo. 


La actividad desarrollada por las células ciliadas parece altamente coordinada. En efecto, 
Assheton (1897) ya había descripto que un embrión de anfibio se puede desplazar por efecto de la 
actividad de sus cilias 1mm/4-7 segundos, siendo el movimiento en dirección cefálica. Es de destacar 
que el DO (que interfiere con la polimerización de al tubulina) inhibe el desplazamiento del em- 
brión por la acción de sus cilias (Herkovits y col., 1975) pareciendo en consecuencia que la actividad 
ciliar estaría vinculada a cambios en el estado de la tubulina. 


En los estadios avanzados del desarrollo embrionario, cuando las células ciliadas involucionan 
se observa que las cilias disminuyen gradualmente su longitud hasta desaparecer. Esta observación 
coincide con la realizada por Kessell y col., (1974) en otro anfibio y sugiere que la desaparición de 
las células ciliadas estaría vinculada a un proceso de reabsorción de la cilias. 


Con respecto a la superficie apical de las células glandulares es de destacar el desarrollo de las 
microvellosidades que siendo numerosas durante el estadio de huevo y blástula (Herkovits, 1976a), 
prácticamente desaparecen dutante el estadio gastrula (Herkovits, 1977b) volviendo a aparecer a par- 
tir del estadio de néurula; ya en estadio 17 las microvellosidades forman una extensa red de pliegues 
sobre la superficie de las células gladulares. Los glóbulos presentes sobre la superficie de estas célu- 
las parecen corresponder a mucus (material secretado por estas células), mientras que los pequeños 
orificios observados sobre la superficie apical podrían corresponder a los lugares donde acontece el 
proceso de secresión. 

Es de destacar que el desarrollo de las células ciliadas y glandulares no parece afectarse en los 
embriones que fueron sometidos a una deformación experimental muy grave durante parte de su 
desarrollo (Herkovits, 1977a, 1978). En efecto, de acuerdo con los resultados obtenidos, dentro de 
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límites muy amplios la forma de las células “per se” no incide en las características específicas de 
la diferenciación celular (v.g. la formación y el mantenimiento de las cilias). 
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FIGURAS 


Figura 1. Células ectodérmicas pertenecientes a la región céfalo-lateral de un embrión de estadío 13,5. (C) célula 
ciliada. 


Figura 2. Células ectodermicas ciliadas y glandulares pertenecientes a la región dorsal de un embrión de estadío 
17,5. (C) células ciliadas agrupadas. 


Figura 3. Ectodermo de un embrión de estadío 21. (C) cilias (G) células glandulares. 
Figura 4. (U) unión apico-lateral entre las células ectodérmicas (C) cilias (G) granulos de secresión. 


Figura 5. Microfotografía obtenida con el microscopio electrónico de transmisión. (U) complejo de unión apico- 
lateral (D) desmosoma (G) granulo de secresión rodeado por membranas. 


Figura 6. Células ectodérmicas ciliadas pertenecientes a la región caudal (aletas) de un embrión de estadío 24. (C) 
célula con población ciliar abundante; (E) célula con población ciliar escasa; (U) célula que presenta una sola 
cilia pudiéndose observar la red de pliegues sobre la superficie apical. 
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SOBRE LAS LLUVIAS EN LA PAMPA HUMEDA Y OTRAS 
REFERENCIAS 


Por Emilio L. Díaz 
RESUMEN 


Se investigan las diferencias de presión en Sud América y mares vecinos en relación con los ciclos solares (A) 
y (B), así como vinculaciones entre las anomalías de las lluvias en la Pampa Húmeda y las variables meteorológicas 
en diversos lugares, de acuerdo al tipo de ciclo. Se calculó una ecuación de regresión múltiple para expresar las ' 
anomalías de las lluvias anuales en dicha región como función de las anomalías del gradiente barométrico sobre el 
Atlántico cerca de la costa argentina, de la temperatura en el Brasil austral, y de la presión en el borde oriental del - 
anticiclón del Pacífico, siendo los coeficientes de determinación 0,60 para los años de ciclos (A) y 0,51 para los 
(B). Una ecuación del mismo tipo se calculó para Valdivia, en base a las anomalías de la presión en Punta Galera y 
del gradiente Santiago-Punta Arenas, con coeficientes de determinación de 0,72 para los años de ciclos (A) y de 
0,59 para años (B). i 

Asimismo, se computaron los desvíos cuadráticos medios de las lluvias anuales en Sud América, correspon- 
dientes a años (A) y (B). 

El estudio de asociación con los niveles de actividad de las manchas solares durante los ciclos (A) y (B) 
señala buena concordancia inversa entre la temperatura en el Noroeste Argentino y la presión en la región de 
Buenos Aires y en el borde del anticiclón del Pacífico, así como coincidencia entre las precipitaciones, temperatu- 
ras y presión en Orcadas y en Stykkisholm (Islandia). 


ABSTRACT 


Difference in pressure over South America and nearby seas in relation to solar cycles (A) and (B) is inves-. 
tigated. A multiple regression equation was calculated to express annual rain anomalies in the crop and cattle 
raising section of Argentina as a function of departures in the pressure gradient near the Argentine atlantic coast, 
temperature in southern Brazil and pressu1e anomalies at the eastern side of the Pacific anticyclon, being R? = 
0,60 for (A) cycles years and 0,51 for (B) years. A study related to levels of sunspot activity for cycles (A) and (B) 
shows inverted agreement among temperatures in Northwest Argentina and pressure over Buenos Aires region and 
at the edge or the Pacific anticyclon; while coincidence is exhibited by precipitation, temperature and pressure at 
the South Orkneys and at Stykkisholm (Iceland). Standard deviations of annual rain anomalies in South America 
for (A) and (B) years are given. 


En trabajos anteriores (1) y (2), entre otros tópicos se analizaron las lluvias y las temperaturas 
en diversos lugares del mundo en relación con los ciclos undecenales del Sol, clasificados en (A) y 
(B) de acuerdo a la distribución de la polaridad magnética en los grupos de manchas. Se denomina- 
ron ciclos (A) los del tipo similar al comenzado en 1913, y (B) los del tipo análogo al siguiente, ini- 
ciado en 1923. 
En este artículo nos referimos principalmente a la Pampa Húmeda argentina, sin perjuicio de 
algunas referencias a otros lugares. | 


1. Aclaraciones. 
1.1. Incluiremos bajo la denominación de Pampa Húmeda el espacio comprendido por el centro y 
norte de Buenos Aires, centro y sur de Santa Fe, y la provincia de Entre Ríos. 


1.2. Los datos utilizados llegan hasta 1974 o 1975 para las estaciones sudamericanas, según los 
casos, y abarcan precipitaciones, temperatura del aire en superficie y presión atmosférica. 


1.3. Por razones de longitud de registros se tomaron, como datos representativos de la serie preci- 
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pitaciones en la Pampa Húmeda, los anotados en Buenos Aires, Concordia y Córdoba (si bien 
esta última queda fuera de dicha región). 


Si se admite que la serie promedio de las precipitaciones observadas en el conjunto Córdoba, 
Rosario, Concordia, Pehuajó, Junín y Buenos Aires es suficientemente representativa de las 
precipitaciones en la Pampa Húmeda, el criterio señalado en 1.3. se justifica en razón de que el 
coeficiente de determinación (1?) entre ambas series es de 0,83 (lapso 1951-1974). 


La correlación recíproca conjunta de las lluvias anuales en Córdoba, Concordia y Buenos Aires 
(1946-1974) arroja un coeficiente de determinación (1?) de 0,18, lo que indica un grado satis- 
factorio de independencia entre los registros de cada estación. 


Los estudios (1) y () fueron completados con el análisis de la presión atmosférica al nivel del 
mar en relación con los mencionados ciclos (A) y (B), figura 1. 


Se introdujo, también, un análisis del comportamiento de ains variables en función de los 
tramos de actividad de los ciclos solares. 


En el curso de este trabajo se utilizarán las siguientes notaciones, que expresan ies anuales, 
promedio de las estaciones que se citan en cada instancia: 


X,: Anom. Lluvias Pampa Húmeda (Córdoba, Concordia y B. Aires); 
X: Difer. Anom. Presión entre la zona de B. Aires y la región Atlántico Sudoeste (B. Aires, 
Montevideo y B. Blanca menos Trelew y Malvinas); 


X3: Anom. Temperatura en la región austral del Brasil (Río de Janeiro, Sao Paulo y Curitiba); 
X4: Anom. Presión en el borde oriental del anticiclón del Pacífico (Santiago y P. Galera); 
Xs: Anom. Temperatura en el Noroeste argentino (La Quiaca, Salta v Mendozá); 

XÍ: Anom. Presión en el Nororete argentino (Salta, Catamarca y Mendoza); 

X,: Difer. Anom. Temperatura entre la Mesopotamia y la Patagonia (Goya menos Sarmiento); 
Xg: Anom. Presión en la zona de B. Aires (B. Aires, Montevideo y B. Blanca); 

Xy: Anom. Presión en P. Galera; 

X10: Difer. Anom. Presión entre Santiago y P. Arenas; y 

X1,1: Anom. Lluvias en Valdivia. 


2. Investigaciones en general 


” 
Le 


Al: 


En el órden anual. 

Los ciclos (A), ver (1) y (2), se presentan un 4 % mas lluviosos que los (B) en el centro y 
oriente argentinos. Asimismo, los años de mayor actividad solar (años 2* a 7” de los ciclos) 
son algo más lluviosos (+ 7 % ) en Entre Ríos y en el Brasil austral que los de menor actividad 
y son, también, ligeramente más cálidos en el centro y noroeste de nuestro país. 
Durante los ciclos (A), las temperaturas en el sur del Brasil (X3 ) exceden un poco a las de los 
ciclos (B) (+ 0,1*C) y, en cuanto al gradiente barométrico sobre el Atlántico Sudoeste, próxi- 
mo a la costa argentina, es algo mas reducido en la primera instancia que en la segunda. Ello es 
congruente con circunstancias favorables a una mayor precipitación sobre la Pampa Húmeda 
(X, ) (mayor frecuencia de frente estacionarios y calientes y afluencia de aire más cálido, y tal 
vez menos estable, desde ei sur del Brasil). 

A continuación informaremos sobre los cómputos efectuados en base a valores correspondien- 
tes a años calendarios. Además, en los cálculos de correlación daremos los coeficientes de de- 
terminación (1? ) por considerarlos más representativos de los vínculos entre variables. 
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La correlación entre (X, ) y (X3) (1885-1974) es de signo positivo y arroja, durante los años 
(A) (N: 44 años) un coeficiente de determinación de 0,44, en tanto durante los años (B) (N: 
48) dicho coeficiente baja a 0,28. 

Con respecto a (X, ) también se anotan correspondencias, ahora negativas, aunque mayores en 
los años (B) que en los (A). Las cifras (1896-1974) de los coeficientes de determinación (X, , 
X> ) son de 0,27 para los primeros (Ng: 44) y de 0,15 para los otros (Na : 37). 

El estudio de la carta (figura 1) de isolíneas de las diferencias de presión atmosférica entre ci- 
clos (A) y (B) muestra una cuña negativa que, desde el trópico, se insinúa en el Noroeste 
argentino, en tanto la presión tiende a ser más alta sobre los océanos vecinos a Sud América 
(+ 0,3 mb). Señalemos, colateralmente, que la comparación entre años de mayor y de menor 
actividad presenta también un area depresionaria (—0,4 mb) en el Noroeste (y otra austral, 
5095). | 
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Por otro lado, los coeficientes de determinación (1901-1974) entre (X, ) y (X¿) son de 0,14 
en los años (A) (N: 33) y de 0.06 en los (B) (N: 44), en el sentido de menores lluvias asociadas 
a excesos en la presión. Una conclusión similar y de signo opuesto se obtiene para las tempera- 
turas (X; ) aun cuando de magnitud algo menor. Las cifras sugieren una relación indirecta. 
Análoga idea surge de la indagación (1901-1974) entre lluvias (X, ) y el gradiente térmico Me- 
sopotamia-Patagonia (X, ), con coeficientes de 0,18 para años (A) y de 0,03 para (B). 

En anticiclón del Pacífico, cuyo borde oriental toca la costa central chilena fue, en razón de su 
influencia sobre los empujes de aire fresco y seco que fluyen desde la Patagonia, motivo de 
consideración. Al respecto, también el grado de relación fue distinto según se tratara de años 
(A) o (B). Los coeficientes de determinación (1899-1974) entre (X, ) y (X, ) son de 0,25 en el 


DE 
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primer caso (N: 34) y de 0,04 en el segundo (N; 44), correspondiéndose un exceso de presión 
en el borde del anticiclón con escasez de lluvias en la Pampa Húmeda. 

En cuanto a la variabilidad de las lluvias anuales en Sud América en relación con los ciclos (A) 
y (B), se computaron los desvíos cuadráticos medios de las precipitaciones de los años calen- 
darios (A) (DCM(A) y (B) (DCm(B), en el lapso 1889-1974. En la o 2 se han trazado las 
isolíneas del mais (DCMm(A)/A(DCM(B). 
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En la misma, el centro y norte de Buenos Aires, y Santa Fe y oriente de Córdoba, presentan 
un área de máxima variabilidad de las lluvias durante los años (A), superior en un 40 % a la de 
los años (B). En cambio, en la parte central del sur chileno (y posiblemente en la Cordillera ve- 
cina) la variabilidad es mayor durante los años (B) (aprox. 25 %). 

Si se tiene en cuenta por ejemplo, que las cifras de Buenos Aires, Goya y Córdoba, se fundan 
en 44 años (A) y en 44 años (B) ello, unido a su men induciría a pensar que la diferencia 


- de variabilidad observada es real. 


Con relación a las variables vinculadas a las lluvias en la Pampa Húmeda, también la variabili- 
dad, medida por los desvíos cuadráticos medios, es mayor en los.años (A) que en los (B); para 
(X, ) es 0,852 y 0,803 respectivamente, para (X3) es 0,568 y 0,467, y para (X4 ) es 0,628 v 
0,589. Esto es congruente con lo anterior. 

En el órden mensual i 

Es interesante notar que cuando se consideran en el árden mensual (tilde) las vinculaciones 
entre las lluvias en la Pampa Húmeda (X”, ) y las temperaturas en el sur del Brasil (X” ), las co- 
rrelaciones señaladas en el órden anual vuelven a manifestarse. En la figura 3 se han graficado 
las diferencias entre los promedios mensuales pertenecientes a los ciclos (A) y (B), en ella pue- 
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de verse el paralelismo de los trazados. 
En cuanto a la correlación entre las anomalías (X”, ) y (X”,) correspondientes a los meses 
calendarios, el cuadro que sigue contiene los coeficientes de determinación: 


y d meses de 
Epoca del año do a N (meses) anos (B) N (meses) 
Octubre a marzo 0,48 60 0,18 59 
Abril a septiembre 0,18 60 0,04 59 


Los meses con mayores precipitaciones, esto es los meses más cálidos, son los que presentan 
vínculos más estrechos entre las variables y este es más grande durante los años (A) que duran- 
te los (B). 

El cómputo de la variabilidad de las anomalías mensuales de las precipitaciones (X”, ) y de las 
temperaturas (X”3 ), revela que los desvíos cuadráticos medios (DC*m), tanto para las unas co- 
mo para las otras, son asimismo mayores en los años (A) que en los (B), guardando correspon- 
dencia entre ambas variables. A continuación y calculados en base a un número de meses igual 
al del cuadro anterior, se indican dichos valores: 


Octubre a marzo Abril a septiembre 
meses (A) meses (B) meses (A) meses (B) 
DC'm An.Llv. P. Húmeda (mm) +1:46 +42 2.50 el 
DC'm An.TT z. Río, S. P., Cur. (€ C) + 1,00 + 0,84 + 1,70 + 1,00 


DIFERENCIA ENTRE LOS VALORES MENSUALES MEDIOS 
CICLOS A= CICLOS BB (%C y mm) (h/4974) 


Lluvias P.HUMEDA 
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3. Ecuaciones de regresión múltiple 


sal 


Lluvias anuales en la Pampa Húmeda. d | 

El estudio de las ecuaciones de regresión múltiple, destinadas a expresar las anomalías de las 
lluvias anuales en la Pampa Húmeda (X, ) en función de otras variables meteorológicas y de los 
tipos de ciclos solares, lleva a concluir que los agentes más significativos al respecto son: 


— la diferencia entre las anomalías de la presión en la zona de Buenos Aires y las del Atlántico 
Sudoeste (X, ); 


— las anomalías de la temperatura del aire en la región austral del Brasil (X; ); y 


— las anomalías de la presión en el borde oriental del anticiclón del Pacífico (X, ). 


Como ya se dijo, es posible suponer una vinculación causal entre dichos agentes y las lluvias en 
la Pampa Húmeda. Las anomalías en el gradiente meridional de la presión se ligan a una mayor 
o menor frecuencia de vientos del este sobre la región y por ende del tipo de frentes asociados, 
la temperatura del aire en el sur del Brasil se liga a la calidad del que afluye hacia la Pampa Hú- 
meda, y las anomalías de la presión en el borde oriental del anticiclón del Pacífico pueden co- 
rresponderse con la frecuencia de los empujes de aire seco y fresco. 

El coeficiente de determinación múltiple es mayor (R?: 0,603) durante los años calendarios 
(A) (N: 34 años) que durante los (B) (R?: 0,510 y N; 44). Las mayores lluvias se asocian a 
gradientes barométricos meridionales menos fuertes que lo usual, a temperaturas del aire supe- 
riores a la normal en el sur del Brasil y a defectos en la presión en el borde del anticiclón del 
Pacífico. La recíproca es cierta. 

Si bien el signo de la acción de los mencionados agentes es igual en ambos casos (A y B), la 
incidencia en la expresión de las anomalías de las lluvias es diferente según se trate de años (A) 
o de años (B). Durante los (A) el control más significativo corresponde a las (X, ) temperaturas 
en el Brasil austral (32 %) y al anticiclón (X4 ) del Pacífico (18 %). Durante los (B) contribu- 
yen principalmente el gradiente meridional (X, ) de la presión (25 %) y las (X; ) (con 22 %). 
Cuando se computan las tres variables referidas, los eventos en el Noroeste argentino (tempe- 
ratura y presión) parecen no tener incidencia ya que las correlaciones individuales con las 
lluvias en la Pampa Húmeda (del órden de r? : 0,09 y al parecer indirectas) son absorbidas por 
las anteriores. | 

En control del borde del anticiclón del Pacífico (X4 ) sobre dichas precipitaciones (X, ) parece 
ser real, en particular para los años (A), pues al suprimir esta variable el coeficiente de determi- 
nación de la fórmula de regresión se reduce de 0,603 a 0.478. Durante los años (B) el efecto 
de la supresión es menor, de 0,510 20479. 

La expresión de las fórmulas de regresión (1899-1974) es: 


— años (A): X, (mm) =-62 X, + 176 X, -119 X, R?; 0,603 
— años (B): X (mm) =-—92 X, + 146 X, — 46 X, RESOLO 


Finalmente, señalemos que el desvío cuadrático medio (DCm) de las (X, ), anomalías de las 
lluvias anuales en la Pampa Húmeda,es de + 207 mm en los años (A) y de + 154 mm en los 
(B). La diferencia es significativa. 

Como ya se manifestó, en general las correlaciones entre las distintas variables (y los diferentes 
lugares) parece más estrecha durante los años (A) que durante los (B), no obstante que en ello 
puede tener algún grado de influencia el número de pares de valores en las correlaciones (34 o 
33 para los (A) y 44 para los (B). 


SOBRE LAS LLUVIAS EN LA PAMPA HUMEDA Y OTRAS REFERENCIAS 3) 


La circunstancia citada se manifiesta, entre otros, en los siguiente coeficientes de correlación 
simple (1901-1974): 


años (A) años (B) 
X. con X, (An.TT Noroeste con An.TT Río, S.P., Curit.) + 0,45 019 
X, con X, (An.PP Noroeste con An.Llv.P.Húmeda) — 0,38 — 0,14 
X¿ con X, (An.PP Noroeste con An.PP P. (Galera) + 0,73 + 0,36 
X, con X¿ (An.PP. P. Galera con An.PP. z.Buenos Aires) + 0,60 -+.0,32 
X, 


, con X, (An. Llv. Valdivia con An. Llv. P. Húmeda) + 0,42 07 


La correlación entre las anomalías de la presión en P. Galera (Xy) y las precipitaciones anuales 
en Valdivia (X,1) es igual (— 0,77) sea que se trate de años (A) como de (B). 


3.2. Lluvias anuales en Valdivia (Chile) 


Las anomalías de las lluvias (X,,) en Valdivia pueden expresarse mediante una fórmula de re- 
gresión en la que intervienen las anomalías de la presión en P. Galera (Xy ), lugar próximo a Valdivia, 
y las anomalías del gradiente barométrico entre Santiago y P. Arenas (Xp). 

El vínculo causal puede atribuirse a que los desvíos de la presión se ligan a la mayor o menor 
frecuencia de perturbaciones en la zona (y eventual corrimiento de anticiclón) y a que las variacio- 
nes del gradiente mencionado significan un incremento o disminución del flujo del aire húmedo del 
Pacífico. 

Aquí también el grado de control de las variables señaladas sobre los desvíos de las lluvias en 
Valdivia (X,1) difiere según se trate de años correspondientes a ciclos (A) o a los (B). En el primer 
caso (N : 33) el coeficiente de determinación es de 0,724, en el segundo (N : 44) de 0,594; en ambos 
el agente de mayor incidencia es la anomalía de la presión en P. Galera (62% y 59 %), observándose, 
con relación al gradiente barométrico meridional, que su efecto, si bien no muy grande, es positivo 
(11%) en los años (A) y ligeramente negativo (0,5 %) en los (B). 

Dicha circunstancia se acusa asimismo, aún cuando debilmente, en las correlaciones simples 
entre el referido gradiente y las lluvias en Valdivia. Cabe anotar que durante los ciclos (A) la presión 
en P. Galera es 0,25 mb superior a los (B), en tanto en Santiago es casi igual en ambos ciclos. 

Las fórmulas de regresión son (1901-1974): 


—años (A): 6 O RIC O/704 
—años (B): X11 (mm) == 371 Xy —10X40) R?: 0,594 


los desvíos cuadráticos medios de las lluvias(X,,)son de + 366 mm en los años (A) y de + 477 mm 
en los (B). ¡ 


4. Asociaciones con los tramos de la actividad de los ciclos solares 


Cuando la temperatura, la presión y las precipitaciones se analizan con relación a los tramos de 
los ciclos solares, independizándolas del curso del tiempo y de las variaciones en la duración de los 
ciclos, se manifiestan, al menos en el área sudamericana y mares adyacentes hasta la región antártica, 
con evoluciones en general interrelacionadas y coherentes. | 

A tales efectos, los tramos precitados fueron definidos en base al curso normal del número de 
manchas, considerándoselos en las siguientes instancias de actividad (A o B): 
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— creciente: años O a 4”, notación 2; 
— máximo: años 2? a 7%, notación 4.5; 
— decreciente: años 52 a 10%, notación 7.5; y 


— mínimo, O tránsito al ciclo siguiente: años 8%, 9% y 10% del ciclo en terminación más los 
años O y 1? del ciclo en comienzo, notación 10. 


Como resultado de ello, cada uno de los puntos que definen la marcha de una determinada va- 
riable en función de los tramos de la actividad solar, es el resultado de promediar los registros de los 
años calendarios pertenecientes a un número de ciclos en el lapso que abarcan los tramos referidos, 
es decir una cantidad de registros que segun los casos oscila entre 15 y 30 años calendarios por pun- 
to. 


EVOLUCION EN FUNCION DE LOS TRAMOS DE ACTIVIDAD DE LOS 
[ICLOS SOLARES Ay B 


Fig. 5 
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En la Figura 4 se muestra la relación inversa entre las temperaturas en el noroeste argentino y 
la presión en la zona de Buenos Aires y en el borde del anticición del Pacífico. 

En la Figura 5 aparece la similitud de la marcha de las precipitaciones, temperaturas y presión 
en Stykkisholm (Islandia, 652N) y en Orcadas (6195), en función de los tramos de la actividad solar 
en los ciclos (A) y (B). Dichas estaciones se encuentran en situación geográfica aproximadamente 
simétrica. : 

La exactitud y corrección de tales vinculaciones depende de contar con registros suficiente- 
mente prolongados; en tal sentido cabe mencionar que cuando se fraccionan las anotaciones de luga- 


res que las tienen de larga data (por ejemplo Roma y Milán) a menudo se presentan divergencias 
entre las fracciones. 


Buenos Aires, 20 de junio de 1978 


(1) “Lluvias y temperaturas relacionadas con los ciclos undecenales del Sol clasificados según polaridad de las man- 
chas”. Anales de la Sociedad Científica Argentina, enero-marzo de 1974. 


(2) “Lluvias y temperaturas asociadas con los ciclos undecenales del Soi considerados según la polaridad de las 
manchas”. Anales de la Academia Nacional de Ciencias de Buenos Aires, 1976. 
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SODALITA DEL NOROESTE ARGENTINO. 
IMPLICANCIAS CRONOLOGICAS Y CULTURALES. 


por Jorge FERNANDEZ! y Margarita MENZEL? 


RESUMEN 


En sitios arqueológicos del norte de Jujuy se han presentado últimamente abundantes cuentas de collar tra- 
bajadas en un mineral que, una vez analizado, resultó ser sodalita. En la Argentina no se conocen yacimientos de 
sodalita, pero en cambio existen en cerro Sapo, Cochabamba, Bolivia, que es el único de toda Sudamérica, y dista 
570 km de la frontera con Argentina. Las minas de cerro Sapo fueron explotadas ya en la I y Il etapa de Tiwanaku 
y con mayor seguridad en las siguientes, hecho que se comprueba por el hallazgo de cuencas de ese material en las 
ruinas de Tiwanaku. Al norte de Jujuy y sud boliviano, la sodalita llega seguramente como un producto del comer- 
cio aborígen. Sin embargo, el hecho de que la sodalita argentina se presente asociada a la cerámica tipo San Pedro 
negro pulido, correspondiente a la etapa II de la cultura de San Pedro, fechada en 250 años D.C., planteá la proba- 
bilidad de que ese tráfico no haya sido directo, sino a través de los oasis del desierto del norte de Chile. 


ZUSAMMENFASSUNG 


In mehreren archáologischen Fundorten im Nordwester Argentiniens (Provinz Jujuy) wurden in der letzen 
Zeit Schumckstúcke gefunden unter denen sich auch solche sus einem blauen Mineral befanden, das spótherin als 
Sodalith bestimmt wurde. In Argentinien sind keine Lagerstátten dieses Minerals bekannt, aber 570 Km weiter 
nórdlich, auf bolivianischen Gebiet, befindet sich die Cerro Sapo-Lagerstáte, welche schon in den Etappen 1 und II 
der Tiwanaku Kultur abgebaut wurde. Diese Tatsache geht aus den Funden dieses Minerals in den Ruinen von 
Tiwanaku hervor. Es gelangte dann sicherlich durch Eingeburenen-tauschuhandel nach dem Siiden Boliviens und 
Jujuy. Da die Sodalith-=stúcke in Argentinien zusammen mit der schwarzen, polierten San Pedro-Keramik gefunden 
wurden, die der Etappe II der San Pedro-Kultur (250 Jahre n. Chr.) zugeschrieben werden, lásst sich annehmen, 
dass dieser Handel nicht direkt verlief,: sondern auf dem Umweg ¡ber die Oasen der nordchilenischen Wúste. (M. M.). 


INTRODUCCION 


Rocas semipreciosas (turmalina, sodalita, lapislazuli, turquesa, etc). empleadas por los pueblos 
prehistóricos de la región andina austral en la elaboración de cuentas de collar y otros adornos, cons- 
tituyen un camino perfectamente tentable por la investigación que tenga como meta la reconstruc- 
ción de incidencias socioculturales por ellos atravesadas. Reina bastante confusión, lamentablemente, 
en lo relativo a la identificación de las materias primas utilizadas pues, no habiéndose superado hasta 
ahora la etapa de su apreciación macroscópica, han sido genéricamente englobadas dentro de las ““tur- 
quesas” y de las “malaquitas”, con lo que su valor indicativo potencial se ha perdido. Confabula en 
esto, principalmente, la identificación verdaderamente difícil de estas piedras, cuya coloración carac- 
terizante es sumamente variable en razón misma de su composición mineralógica compleja, tanto co- 
mo por la extrema rareza de algunas de las especies minerales involucradas. Intentamos abrir con es- 
te trabajo una etapa revalorativa del asunto, sobre la base de la aplicación de técnicas físico-químicas 
modernas y también con una interpretación de la dinámica espaciocultural acorde con éstas. Enten- 
demos que su importancia trasciende el marco meramente mineralógico y aún geológico que pudiera 
conferirsele —nada despreciable, por cierto—, y se aproxima más al de las ciencias sociales por el cú- 
mulo de inferencias novedosas que su estudio posibilita establecer en ese campo específico. 

Dentro de las piedras más raras e indudablemente más hermosas, susceptible de alcanzar un 
bello pulimento por abrasión superficial, se encuentra la sodalita. Se trata de un silicato alúmino-só- 
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dico natural, de color azul de ultramar, a veces con tintes violáceos-vinosos muy suaves. Su escasez, 
es tal, que meritúa sea considerada como una verdadera piedra preciosa: apenas se conocen unos 
pocos yacimientos dispersos por los rincones del globo, tales como los de Kangerdluarsuk en Groen- 
landia, Transilvania y Portugal. El hecho de haber hallado hace poco, en el extremo norte de Jujuy 
y regiones limítrofes de Bolivia, cantidades significativas de cuentas de collar de manufactura indí- 
gena trabajadas en una roca macroscópicamente similar a la verdadera sodalita, y atendiendo a la 
importancia no sólo mineralógica, sino también arqueológica y cultural que de un hallazgo tal se 
desprende, nos ha esforzado en procurar una identificación capaz de quedar por encima de cual- 
quier duda, la que pasamos a reseñar de inmediato. 


DIAGNOSIS 


por Margarita Menzel 


Mediante análisis microquímico (ataque con ácido nítrico y posterior evaporación), se deter- 
minó la presencia de los iones sulfato y cloruro. 

El posterior estudio del mineral por difracción de rayos X (método Debye-Scherrer), permi- 
tió identificarlo como sodalita por comparación con la ficha ASTM correspondiente. La presencia 
del ión sulfato se exterioriza en una ligera desviación de los picos hacia valores menores en 0.02 
Angstrom de los valores teóricos. 


INTERPRETACION ARQUEOLOGICA 


Por Jorge FERNANDEZ 


Bodenbender (1899 : 1925), hizo la primera referencia a la existencia de sodalita en la Argen- 
tina, ciertamente dudosa, ya que la procedencia del espécimen en consideración era desconocida. 
Algo más tarde, Boman (1908, II : 625-631), aporta datos mucho más concretos y referidos a ““piéces 
d'enfilage”” descubiertas por él en diferentes lugares del occidente jujeño, de las cuales una mínima 
parte estaban confeccionadas en sodalita. Boman llama la atención sobre la escasez de esta última, 
ya que se encontraban de 50 a 100 cuentas de “turquesa” por cada dos o tres de sodalita y , además, 
sobre su aparición discontinua en las ruinas de poblados prehispánicos del occidente saltojujeño. No 
encontró ningún ejemplar en Puerta de Tastil, y en Tastil mismo, sólo una. Eran muy raras en Moro- 
huasi, así como en el Pucará de Rinconada. Más frecuentes eran en Gólgota, Queta y Sansana. En 
todos los casos citados por Boman, se trata de cuentas muy pequeñas, discoideas, con perforación 
central y sin duda destinadas a formar sartas. Su diámetro era de aproximadamente 0,014 m y la al- 
tura de sólo 0,10 m. En el Noroeste, estas cuentas reciben el nombre de chaquiras y guaycas (3). Co- 
mo dato próximo, diremos que también von Rosen (1924 : fig. 200), hacia 1901, encontró 18 cuen- 
tas discoideas y solo una cilíndrica trabajadas en sodalita, en el valle de Tolomosa, cerca de Tarija, 
Bolivia, aunque sin haber encontrado ningún ejemplar en su recorrido por la Puna jujeña, en un iti- 
nerario coincidente con el de Boman. La exactitud de las determinaciones mineralógicas está fuera 
de dudas, en el primer caso cubierta por la autoridad de Lacroix, y en el segundo por la de Percy 
Quensel —petrógrafos de primera línea ambos—, quienes practicaron estudios microscópicos sobre . 
cortes delgados y análisis químicos de control sobre los ejemplares de Boman y von Rosen, respec- 
tivamente. 

Carecemos de informaciones posteriores provenientes de otras fuentes, pero según las nuestras, 
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donde realmente es abundante la sodalita es en los poblados prehispánicos en ruinas del extremo 
norte jujeño, en un área que abarca Calahoyo, Cienaguillas y Santa Catalina, asi como las cabeceras 
del valle de Caluni en el Calahoyo boliviano (departamento Sur Chichas).No solo es allí abundante, 
sino que el tamaño de las cuentas resulta excepcional, llegando a tener hasta 50 mm de longitud y 
14 mm de diámetro. Como se comprende, el valor estético de estas cuentas está en razón no sola- 
mente de la cantidad de materia prima empleada, sino también en las dificultades técnicas vencidas 
para su ejecución. En el caso de las cuentas cilíndricas muy largas, la perforación longitudinal central 
-—de casi 3 mm de diámetro—, muy difícilmente se habrá efectuado sin correr el riesgo de quebrar la 
pieza. Se trata de perforaciones bicónicas y cilíndricas, que implican el ataque sucesivo desde ambos 
extremos y que aciertan a unirse en el centro. El pulimento de las superficies se ha obtenido por fro- 
tamiento intenso, por lo que el cuerpo cilíndrico de las cuentas ofrece un tenue facetamiento, visible 
únicamente bajo condiciones particulares de iluminación. Es muy difícil acertar con el tipo de abra- 
sivo utilizado, gue seguramente ha tenido que ser una materia pulvurulenta muy dura, del tamaño de 
varios cientos de mallas. 


Ante estos hallazgos de sodalita manufacturada en el norte jujeño, de inmediato se piensa en 
el problema de su procedencia. Sabemos que esta sustancia mineral se relaciona genéticamente con 
la sienita nefelínica, roca que hasta el presente no se conoce en Jujuy ni en las regiones que geográ- 
ficamente le son próximas. Esto significa que no puede tratarse de un mineral local y, menos aún, o 
por consecuencia, puede serlo su manufactura. Con mayor seguridad ha de tratarse, entonces, de un 
producto de intercambio resultante del comercio verificado entre grupos indígenas del pasado. Ha- 
bría que orientar las pesquisas, por lo tanto, hacia centros en los que la sodalita haya gozado en 
algún momento de difusión y preferencia en grado capaz de generar una moda exportable. 


Fischer (1880 : 370), estudió e identificó como sodalita a una roca utilizada en la confección 
de cuentas de. collar y pequeños idolitos, encontradas por los arqueólogos Stibel y Reiss en las rui- 
nas de Tiwanaku (Bolivia). Al año siguiente, esa determinación fué reiterada por Bamberger y Feus- 
sner (1881 : 580). En Tiwanaku había no solamente cuentas de collar, sino también desechos de fa-' 
bricación, pero la sodalita no era indígena de allí, como tampoco de Ancon, lugar este último donde 
se había hallado en pequeña cantidad por uno de los arqueólogos antes nombrados. Era evidente que 
los trozos de mineral había llegado en bruto de Tiwanaku, y que sólo allí habían sido manufactura- 
dos por los artesanos; pero el yacimiento de procedencia permaneció desconocido mucho tiempo. 
Sólo medio siglo más tarde el geólogo Ahlfeld pudo efectuar el redescubrimiento de las viejas minas 
emplazadas en una de las laderas del cerro Sapo, serranía de Palcaya, provincia de Ayopaya, Cocha- 
bamba, Bolivia. El estudio pertinente, altamente demostrativo, fué efectuado juntamente por Ahlfeld 
y Mosebach (1935 : 414). Geológicamente, se trata de un stock de sienita nefelínica implantado en 
pizarras del Paleozoico inferior, cruzado por un filón gegmatítico de ankerita-sodalita de casi 5 m de 
potencia”. Posteriormente, Ahlfeld (in Ahlfeld y Schneider, 1964 : 338), ha publicado una fotogra- 
fía impresionante de las viejas labores, así como un plano muy detallado con su ubicación. Final- 
mente, Ahlfeld y Wegner (1931 : 286-296), comprueban que la sodalita existente en las ruinas de 
Tiwanaku tiene su foco de procedencia en las viejas minas preincaicas del cerro Sapo, aunque lamen- 
tablemente, sin darnos una seguridad absoluta de que la composición química de la sodalita de Ti- 
wanaku y la del cerro Sapo fuese la misma. Del asunto se ocupó Brendler (1934 : 28-31), cuyos aná- 

lisis químicos, salvo ligeras variantes con los de Bamberger, establecen la identidad de ambas sodali- 

s*. Ahlfeld considera que las minas de sodalita fueron explotadas durante la vigencia de la cultura 
de Tiwanaku, olvidadas durante el desarrollo del imperio incaico —los incas no parecen haberse 
interesado por piedras azules de este tipo, u orientaron su comercio hacia otras regiones—, y reex- 
plotadas por los jesuítas hasta su extrañamiento en el siglo XVIII. Además, en su trabajo con Wege- 
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ner, enfatiza acerca de la sugestiva dispersión de la sodalita hasta el norte argentino, sobre la base; 
de los ya citados hallazgos de Boman. 

Ponce Sanginés (in Ponce Sanginés y Mongrovejo Terrazas, 1970 : 64-65), nos ha enseñado la 
estrecha relación existente entre el dominio tecnológico de diferentes especies líticas y las etapas 
evolutivas por las que atravesó la cultura de Tiwanaku. Con referencia a la sodalita, nos dice que su 
color azul “con seguridad atrajo la pupila del hombre altiplanico desde muy temprano”. Establece 
que en los entierros de la etapa aldeana de Tiwanaku (I, fechada radimétricamente hacia el año 237 
antes de Cristo), ya se hace presente la sodalita. Teniendo en cuenta la estrechez del área territorial 
en ese momento cubierta por la cultura tiwanakota (etapa Il, año 43 D.C.), es seguro que sólo pudo 
acceder a ella a través del trueque. Sólo al alcanzar fisonomía urbana, en la época IM (133-374 D.C.), 
pudo Tiwanaku tener dominio sobre mayores y más importantes recursos rétreos —entre ellos, la 
sodalita—, dominio con seguridad consolidado en la etapa IV (374-724 D.C.). 

De modo que si establecemos que la sodalita presente en toda la región andina austral tiene 
un único centro de procedencia en el cerro Sapo (Cochabamba, Bolivia), puede operar perfectamen- 
te como reloj arqueológico, y en tal caso su aparición en un sitio arqueológico determinado puede 
estar significando: 


12) Que corresponde a un momento anterior a la expansión de Tiwanaku (antes del año 350, 
aproximadamente), lapso en que las minas estaban en posesión de una etnía por ahora no determina- 
ble, pero no tiwanakota. 


2?) Que corresponde a un momento de expansión o de dilatación imperialista de Tiwanaku 
(entre el año 350 y 1050 D.C.), y en tal caso haber arribado como producto directo de esa influen- 
cia cultural irradiante, por estar las explotaciones del cerro Sapo bajo su control directo. 


3) Que corresponde a un lapso comprendido entre los siglos XVII y XVII, durante cuyo 
transcurso los sacerdotes jesuítas hacen explotar dichas minas. 


Evidentemente, las opciones enumeradas vertebran en un sistema lábil por demás, y de inme- 
diato se piensa en un elemento taxativo que posibilite un control más rígido del conjunto. Esta posi- 
bilidad ha de relacionarse con otros elementos arqueológicos asociados a la sodalita. Desafortunada- 
mente, no siempre se presenta ésta asociada a otros restos culturales, siendo lo común que se mani- 
fieste a través de hallazgos superficiales y totalmente desprovistos de vinculación cultural. Además, 
es preciso tener en cuenta la posible reutilización de este material en momentos posteriores, en que 
la materia prima hubiese escaseado por estar cegado el acceso a las minas, sea por guerras u otros 
desacuerdos en el espacio intermedio. Sin embargo, es posible testificar por lo menos un caso en que 
las cuentas de sodalita se han presentado en un entierro dentro de una vasija de cerámica, que era la 
asociación que buscábamos. Ahora bien, la cerámica en cuestión no corresponde a un tipo tiwana- 
kota o tiwanacoide, ni siquiera a alguna de las modalidades de cerámica características de la que- 
brada de Humahuaca o Puna argentina, sino a un tipo cuyo centro de difusión mayor rodea la re- 
gión de oasis de San Pedro de Atacama, El Loa, Chile. Se trata de la cerámica tipo San Pedro negro 
pulido. La aparición de esta cerámica en el norte jujeño y sector sudboliviano limítrofe, así como 
las circunstancias de hallazgo, fueron ya dadas a conocer en dos trabajos (Fernández, MS 1976, MS 
1977), actualmente en curso de publicación. Esta cerámica temprana, cuya edad calendario está 
fijada en 250 años + 150 D.C. (sitio Quitor 6), según Núñez A. (1976 : 74-75), corresponde a la 
fase II de la cultura de San Pedro. En el sitio del norte de Chile indicado (Quitor 6), se presentó 
asociada a componentes de Tiwanaku Clásico y Expansivo (Núñez A., Op. cit, loc, cit, ). Su pre- 


SODALITA DEL NOROESTE ARGENTINO. IMPLICANCIAS CRONOLOGICAS Y CULTURALES 41 


sentación imprevista en el norte de Jujuy, plantea ahora un interrogante sumamente importante por 
las nuevas perspectivas que posibilita avisorar a la investigación. Pero, más sugestivo aún, es que esta 
corporización de la Cultura de San Pedro en su fase Il, se materialice en asociación con un elemento 
como las cuentas de sodalita, tan características en Tiwanaku mismo, y que por nosotros es consi- 
derado como su centro difusor. Desconocemos si en San Pedro de Atacama se ha encontrado soda- 
lita manufacturada; pero, en la situación actual, la asociación apuntada daría pie a inferir que los 
nexos entre el sector del Noroeste argentino aquí considerado y el centro cultural de Tiwanaku 
—fuera de toda duda existentes—, no se han verificado en forma directa, sino a través de un vigoroso 
flujo cultural cuyo foco de irradiación se encontraría en la zona de oasis de San Pedro de Atacama, 
norte de Chile. 
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NOTAS 


(1) Investigader del Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Instituto Nacional de 
Antropología. 


(2) Laboratorio de Mineralogía y Petrografía, Centro Para la Investigación de las Industrias Minerales, INTI (Insti- 
tuto Nacional de Tecnología Industrial). 


(3) Guayca es voz difundida en la región andina. Chaquira, según Zorobabel Rodríguez (1875 : 149), parece ser de 
orígen antillano. Ambas sirven para designar las cuentas y abalorios (granos de alfójar). En la región araucana 
de Argentina (Neuquén y sur de Mendoza), chaquira difunde juntamente con llancatu y llanca; varios Caci- 
ques famosos de los siglos XVI! y XIX la tuvieron por nombre (Llanquitur, Yanquetruz). Tomando en cuenta 
el pequeñísimo tamaño de alguna de estas cuentas, debe pensarse que más bien que a formar collares, lo esta- 
ban a constituir adornos complementarios del peinado y aún de la vestimenta; sólo secundariamente estas sar- 
tas de chaquiras se han utilizado rodeando las muñecas y los tobillos, sectores en los que inmediatamente tien- 
de a emplazarlas nuestro criterio “occidental”. Empero, existen pruebas de su uso alrededor del cuello, Esta 
tradición es sin duda muy antigua en Sudamérica, y ella por sí basta a explicar la increíble aceptación que pos- 
teriormente tuvieron las cuentas de vidrios multicolores, tanto como las venecianas (perlas agri), introducidas 
por los europeos para su tráfico con los indígenas. Igualmente es falsa la clasificación de “elementos de ador- 
no” que de estos objetos tradicionalmente hemos hecho, pues para el aborígen su uso envolvía aspectos mági- - 
co-religiosos para nosotros recónditos. Fundamentalmente, eran objetos sagrados, objetos de ofrenda destina- 
dos no sólo al culto de los muertos, sino también a las maras de la Naturaleza, con las que pr niDan un en- 
tendimiento amistoso a través de estas dádivas. : 


(4) El grupo de los feldespatoides está integrado por la sodalita (Si, O,, Al, Cl, Na, por la noseana (Si, O,, Al, 
Na, J So,) y por la haúyna (Si, Oz, Al, (Na Ca), as *(S0,), y 2) (Klockmann y Ramdohx, 1961 : 642-643). 
Noseana y Haiyna no se encuentran en las rocas profundas, sino que su vinculación genética es con las rocas 
volcánicas jóvenes. La sodalita, en cambio, acompaña o reemplaza a la nefelina en las sienitas. En regiones pró- 
ximas de la Argentina, existe una sola región en la que la sodalita tendría posibilidades de presentarse. Ello es 
en la cordillera Oriental, departamento de Santa Victoria, provincia de Salta. Roberto Beder (1928 : 16-17), 
hace alusión a rocas alcalinas de esa área, y en su mapa geológico ubica dos tifones de sienita rica en sodio. 
Como productos de diferenciación del magma, existen allí filones de rocas variadas (bostonita cuarcífera, bos- 
tonita, basanita, limburgita, etc). Se trata de un conjunto de rocas de veta de extraordinario interés científico. 
En su ya citado mapa, Beder ha delineado un tercer cuerpo intrusivo de edad supuestamente terciaria, en la 
región del río Lipeon, en contacto con areniscas devónicas (área de Baritú, cubierta por selvas de la formación 
tucumano-oranense), aunque sin hacer en el texto ninguna alusión a él. Desconocemos investigaciones poste- 
riores en ese sector; pero desde ya, la asociación de rocas existente en esta última región, y la que sigue inme- 
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diatamente al norte, en territorio boliviano, se presenta como la más favorable para una eventual presentación 
de la sodalita, hecho aún exclusivamente conjeturable, 


(5) Análisis de la sodalita boliviana, según Wolgang Brendler, 1934: 30, 


1 n o IV 


Na, O 24.54 24.53 24.51 22.93 
K,O 1.01 1.13 1.16 0.74 
CaO : - - | 0.46 
Fe, O, -1.06 0.70 0.81 0.85 
Al, O, 31.08 Sa 31.14 30.96 
SiO, 36.63 36.72 36.70 37.96 
Cl 7.20 709 72 5.34 
H, O (110 2C) 0.31 0,23 0.19 1.10 
0 = 0 1.63 1.63 1.63 1.21 


100.20 100.07 100.09 992S 


El análisis de la columna IV, corresponde al practicado por E. Bamberger en 1881. Brendler dice que omitiendo 
el óxido de hierro, estos análisis posibilitan para la sodalita la fórmula 3 Na Al Si O, Na Cl, acorde con los re- 
sultados obtenidos por Linus Pauling (1930)4, trabajando con rayos X. 
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silla de Correo 16. La Plata 1900, Argentina. 


RESUMEN 
Se extienden métodos anteriormente propuestos para el cálculo aproximado de funciones de onda y sus 
autovalores asociados. Se utiliza a tales efectos el teorema no-diagonal del hipervirial y el procedimiento de esca- 
lado óptimo que conduce a la satisfacción del teorema diagonal del hipervirial Empleando como modelo el oscila- 
dor armónico unidimensional, se comparan los resultados obtenidos con los presentados anteriormente y se obser- 
va una sensible mejora en los mismos, 


ABSTRACT 


Previous methods for obtaining approximate wave-functions and eigenvalues associated are extended. The 
non-diagonal hypervirial theoren and the optimum scaling procedure for the fulfilment of the diagonal hypervirial 
theorem are used. The method is applied to the unidimensional harmonic oscillator, the results are compared with 
the previous methods and they show to be in a better agreement to the exact ones. 


1. INTRODUCCION 


Si A y B son autofunciones del operador hamiltoniano H con autovalores E, y Ep respectiva- 
mente, y si W es un operador lineal, acotado en el espacio de Hilbert e independiente del tiempo, en- 
tonces la ecuación de movimiento de Heisenberg en la representación energética establece que[? ] 


(A|[H,W]|B) =(E, — Ep) € A[W|B) 0) 


La ecuación (1) es conocida como el teorema no-diagonal del hipervirial y al operador W se lo 


denomina operador hipervirial. Cuando A y B son degeneradas o idénticas, la ecuación (1) toma la 
forma más sencilla. 


(A]|[H, W]|B> =0 , (2) 


que constituye la expresión matemática del denominado teorema diagonal del hipervirial para el 
operador W [2] y el mismo establece que los valores medios de los operadores que no dependen del 
tiempo no varían con el tiempo en los estados estacionarios. En el caso. particular en que W viene 
dado por 


W=%q:Pi | (3) 


donde la suma se extiende sobre los 3n grados de libertad para el sistema compuesto por n partículas 
y donde q y p denotan las coordenadas y momentos asociados, la ecuación (2) toma la forma 
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1) 
AN 


T y V son. respectivamente, los operadores energía cinética y energía potencial, y se hallan vincula- 
dos al operador hamiltoniano H por la relación 


H=T+V | (5) 


La ecuación (4) constituye la expresión matemática del denominado teorema mecánicocuántico del 
virial [2] y en el caso en que la energía potencial es una función homogénea de grado n, se llega a la 
expresión 


AS (6) 


En esta forma, el teorema del virial resulta de gran utilidad para estudiar las energías electrónicas de 
átomos y moléculas[3 ]. 

S0(M)-=041........- , Y) es una función de onda aproximada, entonces multiplicando 
cada vector posición r¡ por un factor de escala apropiado s? , se consigue que la función aproximada 
Q (r) satisfaga el teorema del virial. La fórmula para el factor de escala apropiado viene dada por: 


s=[n(V0)) /2 CTO) o 0) 


donde n es el grado de homogeneidad de la función potencial V, y ( v(1) ) y ( T(1) ) denotan 
los valores medios de las energías potencial y cinética respectivamente, calculados mediante la fun- 
ción aproximada Q (r) sin escalar. La condición de satisfacción del teorema del virial constituye una 
condición necesaria para toda función de onda aproximada, ya que como lo expresa la ecuación (4), 
toda función de onda exacta lo cumple. En varios trabajos se ha estudiado la aplicación del teorema 
diagonal del hipervirial y su relación con otros teoremas usuales en cálculos y desarrollos teóricos de 
la Química Cuántica [ 5 — 17]. En contraposición a ésto, el teorema no-dianonal del hipervirial no 
ha recibido hasta el presente igual consideración y solamente en un número limitado de casos se lo ha 
aplicado para estudiar problemas de interés químico [18 — 22]. 

Basándose en la propiedad de que para una dada autofunción B y su autovalor asociado Ep, A 
es autofunción de H si y solo si la ecuación (1) es satisfecha para un conjunto suficientemente gran- 
de de operadores hiperviriales W, conn=0, 1,...., N, Coulson y Nash han propuesto el empleo 
del teorema no-diagonal del hipervirial como un medio de cálculo para determinar autofunciones y 
autovalores asociados [?1 ]. 

El método consiste en partir del conocimiento de la función de onda exacta correspondiente 
al estado fundamental y su autovalor asociado, tomar una función de prueba que dependa lineal- 
mente de varios parámetros y elegir arbitrariamente un conjunto de operadores hiperviriales en un 
número igual al de parámetros lineales. Entonces los coeficientes lineales y las diferencias de ener- 
gía E, — E, se calculan por emdio de la ecuación (1) (una ecuación para cada operador hipervirial) 
y empleando la condición de normalización para la función de onda aproximada. Este procedimien- 
to se aplicó al modelo del oscilador armónico unidimensional y los resultados obtenidos resultaron 
bastante alejados de los valores exactos y además, marcadamente dependientes de la forma de la 
función de prueba elegida. En un trabajo reciente, Gopinathan demostró que tales resultados se de- 
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bían fundamentalmente al empleo de un factor de escala incorrecto y propuso un procedimiento 
para el escalado óptimo de las funciones de onda obtenidas por Coulson y Nash[?? ]. Mediante esta 
corrección se consiguen mejorar sensiblemente los resultados, pero entonces las autofunciones obte- 
nidas no satisfacen la ecuación (1), como puede verificarse facilmente por sustitución directa. Se lle- 
ga así a la situación de que la mejora obtenida por el escalado óptimo se logra a expensas del cum- 
plimiento del teorema no-diagonal del hipervirial. En la presente se muestra como introducir el 
factor de escala apropiado sin que se deje de satisfacer el teorema no-diagonal del hipervirial y obte- 
ner funciones de onda aproximadas en tal contexto. 


2. APLICACION DEL METODO AL SEGUNDO ESTADO EXCITADO DEL OSCILADOR 
ARMONICO UNIDIMENSIONAL Y DISCUSION DE LOS RESULTADOS 


El procedimiento propuesto consiste en tomar de funciones de onda escaladas optimamente y 
determinar los coeficientes lineales y autovalores correspondientes por medio de las ecuaciones (1) y 
la condición de normalización para tales funciones. A los efectos de realizar una comparación direc- 
ta con los resultados de Coulson y Nash, y de Gopinathan, se utiliza el método propuesto al mismo 
modelo utilizado por tales autores. Empleando igual nomenclatura que Coulson y Nash, la función 
de onda exacta correspondiente al estado fundamental del oscilador armónico unidimensional y su 
autovalor asociado son 


BSexp (2/2) E, =1/2 (8) 
mientras que las funciones de prueba tienen la forma 


je 
A=exp (+kx%) 2 a (9) 
k es un parámetro que diferencia a las distintas funciones y p es un indicador de paridad (1 para es- 
tados impares y O para estados pares). Como aplicamos el procedimiento al segundo estado excitado 
del oscilador armónico unidimensional (estado par), p es igual a cero en nuestro caso. Los operado- 
res hiperviriales elegidos son de la forma 


Wa =P (10) 


Bajo la forma escalada, las funciones de onda de prueba (9) para el segundo estado excitado. toman 
la expresión 


A(axk)=a1/? (a+ boa? x?) exp (ka? x?) (11) 


donde a y b son los coeficientes lineales a ser determinados y a es el factor de escala. En la Figura 1 
se presentan los gráficos de las funciones de onda obtenidas por el procedimiento de Gopinathan y 
las calculadas por el método propuesto en la presente comunicación, para diferentes elecciones ini- 
ciales del parámetro k. En la Tabla I se dan los valores medios de T, V y E, junto con el solapamien- 
to de las funciones aproximadas con la correspondiente función de onda exacta (k = 0.5). Como 
puede observarse, los resultados obtenidos mejoran sensiblemente los calculados por Gopinathan. El 
mismo tipo de resultado se observa para otros estados excitados, de modo que no son presentados 
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aquí. La mejora obtenida por Gopinathan respecto de los resultados de Coulson y Nash se debe 
esencialmente a la introducción de un nuevo parámetro en las funciones de prueba (el factor de es- 
cala), lo cual, naturalmente, les otorga una mayor flexibilidad a las mismas. Las marcadas discrepan- 
cias en los resultados de Coulson y Nash respecto de los valores exactos, evidentemente están vin- 
culadas al empleo de un factor de escala constante para el término x? en el exponencial (k), mien- 
tras que el procedimiento del escalado óptimo introduce un factor de escala prao en el expo- 
nencial que multiplica a k (a?). 

Como conclusión del presente estudio podemos afirmar que los teoremas hiperviriales ofrecen 
un medio apropiado para obtener soluciones aproximadas a la ecuación de Schródinger, en tanto se 
introduzcan en las funciones de prueba los parámetros necesarios. En aquellos casos en que el con- 
texto de trabajo sea el teorema diagonal del hipervirial, tales requerimientos pueden ser satisfechos 
apropiadamente calculando el factor de escala óptimo, lo cual conduce a que las funciones de prue- 
ba los parámetros necesarios. En aquellos casos en que el contexto de trabajo sea el teorema diago- 
nal del hipervirial, tales requerimientos pueden ser satisfechos apropiadamente calculando el factor 
de escala óptimo, lo cual conduce a que las funciones de prueba cumplan el teorema del virial. Esto 
se puede llevar a cabo mediante la metodología presentada por Gopinathan. Cuando el contexto de 
trabajo es el teorema no-diagonal del hipervirial, una posibilidad para obtener buenas funciones de 
prueba es emplear el procedimiento presentado en este trabajo, el cual aparece como muy conve- 
niente para su aplicación a otros tipos de problemas, en vista de la calidad de los resultados obteni- 
dos para el modelo del oscilador armónico unidimensional, que son superiores al del escalado ópt:- 
mo. Los teoremas hiperviriales se constituyen así en una alternativa válida con respecto a los clá- 
sicos métodos variacionales y perturbacionales. Actualmente se está trabajando en nuestro laborato- 
rio en la aplicación del método aquí presentado a otros modelos de interés químico, así como en el 
estudio de la extensión y generalización del mismo. 


Función de Gopinathan!221 Este a 
NE E | de | Solapamiento| i Eos l [Solapamiento 
» [| m | | ade 


0,9749 
0,9963 
1,0000 


0,9984 
0,9952 
0,9918 


Tabla 1. Valores medios energéticos y solapamiento con la función de onda exacta del segundo estado excita- 
do del oscilador armónico unidimensional. 
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Figura |. Funciones de prueba hiperviriales para el segundo estado excitado del oscilador armónico 
unidimensional de la forma dada por la ecuación (11). Los gráficos a línea punteada corres- 
ponden a las funciones de Gopinathan (Ref. 22) y los gráficos a línea llena corresponden a 
funciones determinadas en este trabajo. 
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